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Abstract. We use super-paramagnetic spherical particles which are arranged in 
a two-dimensional monolayer at a water/air interface to investigate the crystal to 
liquid phase transition. According to the KTHNY theory a crystal melts in thermal 
equilibrium by two continuous phase transitions into the isotropic liquid state with an 
intermediate phase, commonly known as hexatic phase. We verify the significance of 
several criteria based on dynamical and structural properties to identify the crystal 
- hexatic and hexatic - isotropic liquid phase transition for the same experimental 
data of the given setup. Those criteria are the bond orientational correlation function, 
the Larson-Grier criterion, 2D dynamic Lindemann parameter, the bond-orientational 
susceptibility, the 2D Hansen- Verlet rule, the Lowen-Palberg-Simon criterion as well 
as a criterion based on the shape factor of Voronoi cells and Minkowski functionals. 
For our system with long range repulsion, the bond order correlation function and 
bond order susceptibility works best to identify the hexatic - isotropic liquid transition 
and the 2D dynamic Lindemann parameter identifies unambiguously the hexatic - 
crystalline transition. 

PACS numbers: 64.60. Q-, 81.10.-h, 82.70.Dd 
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1. Introduction 

While the liquid to crystal transition in three dimensional systems is usually a first 
order transition the situation in two dimensional systems is found to be more complex. 
While grain-boundary induced melting [TJ [2] or condensation of geometrical defects 
[3j H] suggest a first order phase transition, the theory of Kosterlitz, Thouless, Halperin, 
Nelson, and Young [5], [6], [TJ, [8] predicted a melting process via two continuous phase 
transitions with an intermediate phase. The intermediate phase appears due to the fact 
that the translational and orientational symmetries are broken at different temperatures. 
The first phase transition at temperature T m is associated with destroying the discrete 
translational symmetry. The intermediate phase is named hexatic based on the 
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remaining sixfold quasi long range orientational order. If the orientational symmetry 
is destroyed to a short range one at temperature Tj > T m a second phase transition 
from the hexatic to the isotropic liquid occurs. According to the KTHNY theory the 
different symmetries are affiliated with the occurrence of different topological defects, 
namely dislocations and disclinations. 

The first simulations of small systems of hard discs showed a single first order 
transition [9]. Increasing system size a first order transition with short correlation length 
was ruled out but data were compatible with a weak first-order transition as well as a 
continuous scenario. Since a phase coexistence was reported in an equilibrated ensemble 
we would like to argue in favor of a weak first-order scenario. Simulating larger systems 
with up to 4 million particles with the same computer code it was shown that a Van-der- 
Waals loop, which is usually interpreted as a first order criterion, weakens with increasing 
system size [10J. But especially in small systems the existence of a Van-der-Waals loop 
can not be taken as solely criterion since the size of the loop strongly depends on the 
boundary conditions and even systems wich are known to have continuous transitions 
show a Van-der-Waals loop [11]. Simulations of dipolar particles were consistent with 
KTHNY theory [12] . Wether a system melts via a first order or via KTHNY theory may 
depend on the core energy of dislocation [13] a quantity which one can calculate a priory 
only in the dilute limit of dislocations where renormalization effects do not appear [13] . 
Binder et al. pointed out that the KTHNY scenario may easily be preempted by a first 
order transition [15]. Recent large scale simulations of hard core particles reported a 
continuous transition between crystal and hexatic phase but since a phase coexistence 
of hexatic and isotropic liquid was found, this transition is reported to be first order 

HE]. 

In experimental systems, the existence of the intermediate hexatic phase is well 
established pU [HI H3 EH ED 1221 1221 123] but the nature of the transition is debated as 
well. Indications of first order transitions are reported in a colloidal system with screened 
coulomb interactions [17] and di-block copolymer systems [211 [22]- Like in simulations, 
a phase coexistence is usually interpreted as fist order signature. KTHNY theory is a 
melting theory starting from large single crystals but of cause a thermodynamic phase 
should be independent of the history of the matter and cooling and heating cycles 
should serve the same results. Wang et al. reported to find poly-crystalline domains 
during cooling in a system of diameter-tunable microgel spheres at finite cooling rates 
[25]. Indeed, cooling rates have to be small such that critical fluctuations of continuous 
order phase transition can switch the symmetry globally. In the present system of 
particles with dipolar interaction we found (within the given resolution of temperature) 
both transitions to be continuous during melting and freezing. If the cooling rate is 
very slow (keeping the system always in quasi-thermal equilibrium) we do observe large 
single crystal domains in the field of view [23] implying that KTHNY theory also holds 
for freezing. But when cooling the system rapidly from the isotropic liquid to the 
crystalline state we find a poly-crystalline sample without a signature of a hexatic phase 
during crystallization [26]. But even if growing crystalline domains are found in a liquid 
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environment shortly after a temperature quench one should not interpret them as liquid 

- crystal coexistence - simply because the system is far out of equilibrium. In the same 
sense one should be careful taking poly-crystallinity solely as a signature of first oder 
transition if the system is cooled at a finite rate. 

The manuscript is organized as follows: after a short introduction about long range 
order of crystals in two dimensions the melting theory developed by Kosterlitz, Thouless, 
Halperin, Nelson, and Young is summarized in section |2j In section [3] we introduce our 
experimental setup and how we realise a colloidal monolayer. The following sections H] 

- [11] introduce several quantities to identify different thermodynamic phases in 2D and 
discuss the results for our colloidal system with dipolar particle interaction. Finally we 
summarise the advantages and disadvantages of the measures in the conclusion. 



2. Crystals in two dimension and KTHNY 

Since the work of Peierls [27] 128] and Mermin [29] 130] it is known that strictly speaking 
no crystals exist in systems with dimension D < 3. In general the significance of 
fluctuations is increased if the dimension of a system is decreased. Crystal lattices with 
dimension D < 3 are thermally unstable due to long-wavelength phonon modes. As a 
consequence long-range translational order does not exist. In case of D = 2 Mermin 
showed that the displacement autocorrelation function 

([u(r) - u(r')] 2 ) ~ln|r-r'| |r - r'| -»> oo (1) 

diverges logarithmically in the crystalline phase. The slow logarithmic divergence in 2D 
leads to crystals which posses only a quasi long-range translational order. On the other 
hand the local crystalline orientation is preserved and a long-range bond orientational 
order exists. The absence of long-range translational order also affects the shape of the 
structure factor S(q). In 3D the structure factor is characterized by a number of delta 
functions 

S(q) ~ 8(q - G) (2) 

at the reciprocal lattice vectors G reflecting a not diverging displacement u(r). In 2D 
the delta functions are replaced by a set of power-law singularities 

S(q)Hq-Gr 2+rte ( T ) (3) 

where 

k B T\G\H3, R + X H ) 
,Gl ; 47r/ lH (2^ + A R ) 1 ' 

depends on the Lame coefficients and \r. 

As mentioned before, a 2D crystal is characterized by a long-range orientational 
and a quasi long-range translation order. The KTHNY theory describe the melting of a 
hexagonal crystal by the appearance of thermally induced topological defects. The phase 
transition crystal - hexatic takes place at temperature T m when thermally generated 
bounded pairs of dislocations which spontaneously appear in the crystal phase dissociate 
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into single dislocations. While dislocations destroy the quasi long-range translational 
order a quasi long-range orientational order is preserved. The second phase transition 
hexatic - isotropic liquid at Tj > T rn occurs when dislocations are separated into free 
disclinations. Disclinations destroy the remaining quasi long-range orientational order so 
that in the isotropic liquid both the translational and orientational order is short-range. 

By identifying the kind of the translational and orientational order at a given 
temperature T the state of the system can be defined and thus the phase transition 
points. The nature of the order is marked in a different behaviour of the density density 
correlation function go and bond orientational correlation function g§. The density 
density correlation function is given by 

9c{r = \r k - rj|) = (p*(r fe )p(r,)) (5) 

where p(r^) = e tGvk is the translational order parameter of particle k located at position 
rfc. Practically the translational correlation function is rarely used to identify the 
melting temperature T m . This is due to the fact, that in a system with Mermin- Wagner 
fluctuations being present, the reciprocal lattice vector G is not easily determined due 
to the power law singularities of the structure factor. Especially in large systems it is 
numerically difficult to extract G unambiguously. 

The bond orientational correlation function is defined by 

9e(r = \r k - r { |) = (V>g(r fc )^ 6 (r,)) (6) 

with the sixfold bond orientational order parameter 

W = -E^ (7) 

where n\ is the number of nearest neighbors of particle k and Q k \ is the angle between the 
bond of particles k and I and an arbitrary but fixed reference axis. Here the () brackets 
correspond to an ensemble average. The long-range orientational order in a crystal 
is expressed in a long-range bond orientational correlation function lim^oo g§{r) ^ 
whereas an algebraic decay of the density density correlation function gc(r) ~ r~ VG ^ 
reflects the quasi long-range translational order. The hexatic phase is characterized by 
an algebraic decay of g§{r) ~ r~ r>6 ^ (q Uas i i on g rari ge orientational order) with an 
exponent 

18k B T 

Ve{T) = —^r- (8) 

depending on Frank's constant Fa- The short range translational order on the 
other hand leads to an exponential decay of gc( r ) ~ e ~ r ^ G ^ with £g(D being 
the translational correlation length. In the liquid regime of the phase diagram the 
orientational order is short range too, and the correlation function decays as g§{r) ~ 
e -r/£,e(T) w h ere £ 6 (T) is the orientational correlation length. 

In addition to the predictions of the KTHNY theory different criteria have been 
proposed to identify the phase transitions. Whereas Wang et al. [25] tested various 
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Figure 1 . Cutout of the image of the CCD-camera of a monolayer of colloidal particles. 
For high magnetic fields particles arrange in a crystal (left) and for low magnetic field 
particles form a fluid (right). Images are a quarter in size of the full field of view 
(about 580 x 430 /im 2 ). The hexatic phase is hardly distinguishable from the isotropic 
phase solely by eye and ist not shown here. 



2D freezing criteria in poly- crystalline samples of microgel particles we do so in mono- 
crystalline samples of dipolar particles. We verify the adaptability of the bond order 
correlation function, the Larson-Grier criterion, the Lindemann parameter, bond- 
orientational susceptibility, 2D Hansen Verlet rule, 2D Lowen-Palberg-Simon criterion, 
shape factors and Minkowski functionals on the melting transitions. 



3. Experimental System 

The experimental system is described in detail in [3T| . It consists of a 2D colloidal 
monolayer of spherical polystyrene spheres with diameter d = 4.5 /im suspended in 
water and sterically stabilized with the surfactant Sodiumdodecylsulfate. Nanoparticles 
of Fe2C>3 are embedded homogeneously in the polystyrene spheres being responsible for 
super-paramagnetic behavior and a relatively large mass density of 1.5 g/cm 3 . Therefore 
particles are confined by gravity at a water/air interface formed by a droplet which is 
suspended by surface tension in a top sealed cylindrical hole (6 mm diameter) of a glass 
plate. An external magnetic field H perpendicular to the water/air interface induces 
a magnetic moment M = %H in each particle causing a repulsive dipole-dipole pair- 
interaction E magn between them. The dimensionless interaction parameter V which is 
given by the ratio of the magnetic versus thermal energy 



Emagn = ftp {.X^f (lX P ) 3/2 x 

knT 4vr knT ^ sys 



oc T~\ (9) 



vB J- n-B -> 

is equivalent to an inverse system temperature. Under the conditions of temporally 
constant ambient temperature T and 2D particle density p the system temperature 
depends only on the magnetic field. As a result the system temperature can be easily 
adjusted by simply changing the strength of the magnetic field H. 
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Figure 2. Orientational correlation functions ge(r) (in units of particle distances a) 
for different interaction parameters T. For T > 60, g§(r) reflects the long-range order of 
a crystal (upper three curves) while in the hexatic phase an algebraic decay is observed 
(two curves in the middle). An exponential decay is observed in the isotropic liquid 
(lower three curves). 

An inhomogeneous distribution of the particles within the sample would induce a 
gradient in system temperature causing a spatial dependence of the phase transition. 
Therefore it is crucial to align the water/air interface absolutely planar and horizontally. 
For this purpose several computer controlled regulation loops have been installed to ad- 
just the interface and keep it temporally constant. A monochrome CCD camera is used 
to observe particles by video microscopy. The field of view (1158 x 865 fim 2 ) contains 
~ 9000 particles whereas the whole system includes ~ 250000 particles. During data 
acquisition the coordinates of the particles in the field of view are determined in situ 
every w 2 sec over a period of 25 min by digital image processing with an accuracy of 
about 50 nm. This way the phase space information is accessible on all relevant length 
and time scales. A crystal is melted by stepwise increasing the system temperature via 
a reduction of the magnetic field. In the range of the phase transitions the interaction 
parameter was changed in small steps with an increment of Ar 0.25. After each 
modification of the interaction parameter the system was equilibrated for about a day. 
Figure [1] shows images of the colloidal monolayer in the crystalline (left) and isotropic 
fluid phase (right). 
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4. Orientational correlation function 

As key-quantity of KTHNY theory, figure [2] shows the bond orientational correlation 
function for temperatures in the crystalline, hexatic and isotropic liquid regime. The 
evaluation of ge(r) includes a time average over several particle configurations in addition 
to the ensemble average. The bond correlation function stays finite for interaction 
parameters T > Tm = 60 in the crystalline phase. An algebraic decay is observed for 
T = 59.6 and T = 59.3 so that the phase transition crystalline - hexatic takes place 
between T = 60 and T = 59.6. In the liquide state (r = 47.6 and T = 39.5) the 
decay of g§(r) is exponential. This behavior is in accordance with the predictions of the 
KTHNY theory and reconfirms previous experimental results [T9l [20j |23] for the given 
system. The hexatic - isotropic transition can be determined by the investigation of the 
goodness of fit statistics of algebraic or exponential decay [23]. In principle one could 
also determine the hexatic - crystalline transition by investigating the slope of g§(r) in 
a log-log plot but it is numerically not very precise to distinguish between a small but 
finite and zero slope decay. From g§{r) one can also extract two diverging quantities: 
The orientational correlation length diverges at Tj and Frank's constant diverges at 
T m . In [23J we fitted both divergencies to extract critical exponents but the transition 
temperatures were not taken as fitting parameter but used as an input. Otherwise fitting 
the transition temperatures to a single sided divergence in a finite field of view always 
overestimates the transition points by a few per cent. 

5. Local bond-order by Larson-Grier 

To get insight into the local symmetry we focus on the magnitude of the local bond 
order parameter 



m§ k is zero for perfect five- or sevenfold neighbored particles and one for perfect sixfold 
ones. It measures how the neighbors of particle k fit locally on a hexagonal lattice. To 
compare the local sixfold symmetry with neighboring particles Larsen and Grier [32] 
investigated the magnitude of the projection of ipe k 



to the mean local orientation field. It takes the second nearest neighbors into account 
and determines how the orientation of particle k fits into the orientation of its neighbor 
particles. Since it is a projection n 6k < m 6fc and n 6fe + m 6k < 2. In [32] an uni-modal 
distribution was found even if real space images showed a dilute liquid (or gas) phase 
and dense crystalline flakes implying an attractive interaction between particles to exist, 
whereas in [17J a bimodal distribution is reported next to the isotropic-hexatic as well 
to the hexatic-crystalline transition. Particles in the mg-n6-plane with m§ + > 1 
(upper right corner) where identified to be crystal-like particles. Figure [3] shows the 
probability distribution for our purely repulsive system in the m6-W6-plane for several 




(10) 
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Figure 3. Probability distribution of the magnitude of the local bond order parameter 
mg versus the magnitude of the projection of ipe to the mean of nearest neighbors tiq. 
The probability distribution changes continuously at both phase transitions and no 
bimodal distribution can be found. 
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Figure 4. The dynamic Lindemann parameter stays finite in the crystal phase and 
diverges in the hexatic and isotropic liquid state. Y increases from top to bottom with 
the same order as in the label. 

temperatures. The upper row are plots in the crystalline phase, the second row shows 
plots of the hexatic phase whereas the two lowest rows are all from fluid phases. The 
absence of a bimodal distribution and the weak dependence of the local bond-order field 
above and below and T m indicate continuous phase transitions. The third row shows 
that the local order in a 2D fluid is predominately hexagonal, even far away from the 
phase transitions. Only at very high temperatures (low interactions strength) most of 
the particles have m 6 + < 1 indicating that the local six-fold order up to the second 
shell is lost (lowest line). Since the dependence of the local bond order on the different 
phases is weak, it does not serve as a sharp criterion for phase transition temperatures. 

6. Lindemann parameter 

The Lindemann parameter is a well-known criterion in 3D to identify the melting point of 
crystalline structures. According to Lindemann [33] the melting of a crystal takes place if 
the thermal energy leads to displacements of atoms in relation to their equilibrium lattice 
sites which are in the range of one-half of the interatomic distance. The Lindemann 
criterion was modified by Gilvarry [31] by considering the root-mean-square amplitude 
of thermal vibrations. He suggests that the melting process is initiated when the fraction 
of the root-mean-square amplitude and the interatomic distance reaches a critical value 
of approximately 0.1. 

The Lindemann criterion in this form is inapplicable in 2D. Due to long wavelength 
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jL{t) = — 2 (13) 



phonon modes the Lindemann parameter diverges in a crystal as well as in a liquid. 
Bedanov et al. [35] introduced a melting criterion for 2D 

7m = <K^±iQ (12) 

cr 

analog to the Lindemann parameter in 3D based on the displacement Uj of particle j with 
respect to its nearest neighbors j + 1 and normalised to the average inter-particle distance 
a. Zahn et al. [12] generalised equation (TT2"j) to a dynamic Lindemann parameter 

([Aui(t) -Au. i+1 (t)] 2 ) 
2a 2 

where Au(t) = u(t) — u(t = 0). The crystal - hexatic phase transition can be determined 
by the long time behaviour of 7z(t). The Lindemann parameter 7^(£ — > oo) diverges in 
the hexatic and isotropic liquid phase whereas in a crystal jjj(t — » oo) stays finite below 
a critical value 7£ = 0.033. 

The Lindemann parameter 7l(£) is shown in figure (jl]) for different interaction 
parameters T. As expected 7z(£) converges in a crystal (r > 60) to a finite value < 7£ 
but diverges in the hexatic phase (T = 59.6 and T = 57.7) and isotropic liquid state 
(r = 53.7, T = 49.6 and T = 39.5). According to the behaviour of the Lindemann 
parameter the transition crystal - hexatic occurs in the range between T = 60.0 and 
T = 59.6. This result is in excellent agreement with the predictions of the KTHNY 
theory and the determination of the phase transition crystal - hexatic obtained with the 
help of the bond orientational correlation function g^. If grain boundaries are visible 
due to finite cooling rates or density gradients in the sample, the dynamic Lindemann 
parameter is not finite for the crystalline state. But for a mono-crystalline sample it acts 
as a very sensitive tool to determine the crystal - hexatic phase transition temperature 



7. Bond-orientational susceptibility 

The hexatic - isotropic liquid phase transition is associated with fluctuations of the 
orientational order parameter The fluctuations can be quantified by the bond- 
orientational susceptibility 

X* = A{(\*l\)-(\* 9 \) 3 ) (14) 

where ^ 6 = l/JVE^iV'e 

(r k ) is the global bond orientational order parameter of the 
iV particles included in a system with size A. The bond-orientational susceptibility 
increases dramatically if the temperature reaches the point of the hexatic - isotrop 
liquid phase transition at Tj [36J. Here an increase of \& is observed independently 
whether the system approaches the transition point from the liquid V — > T~ or from the 
hexatic phase r — > Tf . This behavior of the bond-orientational susceptibility simplifies 
the identification of the hexatic isotropic liquid phase transition in comparison with 
the previously mentioned method of the single sided divergence of the orientational 
correlation length T — > iy. In figure [5] we see a sharp increase of Xe at T = 57.5 ± 0.5. 
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This result coincides with the value for the hexatic isotropic transition obtained in |23j. 
The bond-orientational susceptibility is a very sensitive tool to determine the transition 
temperature since we find a sharp increase from both sides of the peak (unlike e.g. 
the divergence of the orientational correlation length calculated from ge(r) where the 
divergence is single sided, see section H] or [23]). In principle one could use the bond- 
orientational susceptibility as criterion for first order or second order transition (36]. A 
symmetric peak shape is predicted for second order or continuous transition whereas for 
first order transitions the limit of the susceptibility from above and below the transition 
should be different. Due to the limited temperature resolution in our data we do not 
want to overestimate this topic but on a first glance the data seem to be consistent with 
a second order transition. 

8. Structure factor 

The structure factor is another often used physical quantity to identify the freezing 
transitions. It is defined by 




(15) 



where the spatial Fourier transform of the number density p(q) is given by 



iV 




(16) 
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Figure 5. The bond-orientational susceptibility for interaction parameters in the 
isotropic liquid, hexatic and solid state. The maximum of the peak corresponds to the 
isotropic liquid phase transition. 




Figure 6. The structure factor 5(q) in the liquid phase (top) and crystalline phase 
(bottom). The rectangular cross in the center is an artifact due to the finite field of 
view. 



Comparison of 2D melting criteria in a colloidal system 



13 



5? 



40 
35 
30 
25 
20 
15 
10 
5 




T 1 1 1 1 1 1 1 1 1 1— 



r = 


79.2 


r = 


63.4 


r = 


60.4 




58.9 




57.2 




52.6 




45.8 




6.0 




Figure 7. The isotropic structure factor S(q) for different interaction parameters. 
Curves are shifted for reasons of clarity, T decreases from top to bottom with the same 
order as in the label. 



Here () denotes an ensemble average over N particles located at positions r^. Since we 
know the time dependent trajectories of the particles in the field of view we calculate 

the structure factor as function of time 

I j? JV 

S(q,t) = — ^^exp(zqr i (t))exp(-zqr j (t)) (17) 
iV i=ij=i 

and determine 5(q) = nj l J2t=i by a time average over n t > 70 statistically 

independent particle configurations. 

Figure M show the structure factor S(q) in the isotropic liquid and the crystalline 
phase calculated from particle trajectories. Performing an azimuthal average gives the 
classical structure factor S(q). Hansen and Verlet argued that freezing is associated with 
the amplitude S(qo) of the first maximum of the isotropic structure factor. A 3D liquid 
freezes when S(qo) exceeds a characteristic value of 2.85 [37]. The predictions of the 
characteristic value in 2D resulting from simulations vary from S(qo) = 4.4 for particles 
with hard core and coulomb interaction [38] to S(qo) = 5.75 [39J for particles with r~ 12 - 
interaction. Figure (J7j) shows the temperature dependent isotropic structure factor S(q) 
which is obtained by an angular average of the structure factor S(q). The maximum of 
the isotropic structure factor S(q ) increases continuously as the temperature decrease 
(increasing interaction parameter T). Additionally the second maximum splits in two 
peaks which reflects an evolving hexagonal structure. As shown in figure ([8]) S(qo) rises 
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Figure 8. Amplitude of the first maximum of the structure factor S(qo). At the 
freezing point r = 60 the amplitude S(qo) ~ 10 for the system with long range particle 
interaction. Two error-bars are shown in the isotropic and crystalline phase, calculated 
as time average from different time steps at the given temperature. 



slowly in the liquid and hexatic phase followed by a sharp rise after crossing the freezing 
point T = 60. The characteristic value is S(qo) — 10 and thus exceeded the estimated 
value by a factor of ~ 2 for our system with long range dipole-dipole interaction between 
the particles. 

Another criterion about the global order is given by the line shape of the angular 
intensity of the structure factor. According to [4Tj| |4"T] the line shape of the Bragg peaks 
in the solid state are given by a Lorentzian function S(9o) = [(9 — 6) 2 + k 2 }^ 1 where 
#o is the angular position of the maximum of a Bragg peak, k is the angular width 
of the Lorentzian function, and the in-plane angle 9 ranges from 9 — 7r/6to# + 7r/6 
because of the sixfold symmetry. In the hexatic phase a square-root Lorentzian (SRL) 
behavior S(9 ) = [(9 — expected. The line shapes of a Bragg peak in 

the solid T = 80.4 and hexatic F = 58.0 phase are shown in figure M whereas in the 
isotropic liquid state T = 39.5 no angular dependence of the intensity is observed. To 
evaluate the behaviour of the line shapes for different interaction parameters in the solid 
and hexatic phase a fit with a Lorentzian as well as a square-root Lorentzian function 
was executed. The line shape was determined on the basis of the reduced chi-square 
goodness-of-fit statistic \ 2 °f the fits. The ratio Xl/Xsrl °f the Lorentzian and square- 
root Lorentzian reduced chi-square goodness-of-fit statistic is given in the inset of figure 
|9j The line shapes in the solid state are well reproduced by a fit with a Lorentzian 
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Figure 9. The line shape of a Bragg peak in the solid (T = 80.4) and the fit with 
a Lorentzian function (dashed line) and in the hexatic phase (r = 58.0) with a SQL 
function fit (dotted line). In the isotropic liquid (r = 39.5) the intensity shows no 
angular dependence. The curves are shifted for clarity. The inset shows the ratio 
of the goodness-of-fit-statistic as function of system-temperature. Above T m = 60 
Loerentzian function fits better, below T m a SQL function fits better (except the data- 
point at r = 58, see main text). 

function Xl/Xsrl < 1 while in the hexatic phase a square- root Lorentzian function 
fits better Xl/Xsrl > 1- Only in the vicinity of the hexatic to isotropic liquid phase 
transition at Tj = 57.5 this is not the case for the data-point at V — 58. This might 
be due to the fact that below T; the peaks should disappear and this datapoint is to 
close to the isotropic transition to distinguish unambiguously between Lorentzian and 
square-root Lorentzian azimuthal shape. 

9. Lowen-Palberg-Simon criterion 

Lowen, Palberg and Simon jl2] introduced a freezing criterion based on the dynamical 
properties of 3D systems. Their criterion states that a system starts to solidify if 
the ratio of the long-time self-diffusion coefficient Dl and the short-time self-diffusion 
coefficient D reaches a critical value of 0.1. Brownian dynamics simulations of different 
pair potentials lead to critical values in 2D between 0.072 (hard disks) and 0.099 (r~ 12 
potential). In case of a dipolar interaction a critical value of 0.086 was obtained |43j. 
The short-time and long-time self-diffusion coefficients are related to the mean-square 
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displacement Ar 2 = l/NJ2f=i[ r i(t) ~ r i(0)] 2 by the following equations: 

Ar 2 

D Q = lim — (18) 
t^o At v ' 

Ar 2 

D L = lim — — (19) 

The ratio of the short-time and long-time self-diffusion coefficients for different 
interaction parameters are shown in figure [T01 The critical value expected for the 
freezing point from simulation with a dipolar interaction can not be reproduced. At 
the solid hexatic phase transition Dl/Dq ~ 0.03 which is a factor of three smaller 
compared to computer simulation for about 1000 dipolar particles [33]. In a poly- 
crystalline system of soft spheres a threshold of about Dl/D = 0.08 is reported 
[25] supporting the simulations. Since we know that grain boundaries influence the 
2D dynamical Lindemann parameter we expect a difference in the long time diffusion 
coefficient between systems with and without grain boundaries. Additionally the short 
time diffusion coefficient of particles at a water/air interface is larger compared to 
particles in the bulk. This together with the fact that grain boundaries are not visible 
in our field of view might explain the small value of 0.03 in our system. Zippelius 
[33] pointed out, that the dependence of the ratio D L /D as function of temperature 
may be used as criterion for first and second order transition. First oder transitions 
are characterized by a discontinuous jump at T m whereas a smooth change should be 
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Figure 10. The temperature dependent ratios of the short-time versus long-time self- 
diffusion coefficients D^/Dq. At the freezing point D^/Dq ss 0.03 which has to be 
compared to a value of Dl/Do — 0.086 expected from simulation. The inset shows 
the time dependent diffusion constant calculated for different time-windows from the 
mean squared displacement. 
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Figure 11. The fraction of particles assigning to three different classes A,B and C. 
Error-bars (calculated as time average for different time steps at a given temperature) 
are smaller than the symbols. 



found for KTHNY-like behaviour. Figure [TO] shows a continuous variation as function 
system temperature. This again supports the result that our system melts according to 
KTHNY theory. 

10. Shape factor 

Moucka and Nezbeda introduced a shape factor ( to analyze structural changes occurring 
for simulations of a hard disk system in the region of phase transitions [45J . The shape 
factor of particle % is defined as 

* = s| (20) 

where and Si corresponds to the perimeter and area of the Voronoi cell of the particle. 
The shape factor for a regular polygon with n edges is given by CJ^ 9 = n/7rtan(7r/n). 
They observed that the distribution of the shape factor P(() becomes bimodal near the 
freezing point. The distribution in the liquid state is broad and the maximum is located 
at relative high (^-values whereas in the solid a sharp distribution with a maximum 
near the value of a regular hexagon ( = 1.03 is given. This behaviour is explained by 
the evolution from different type of cells which are also distorted in the liquid to more 
regular hexagonal cells in the solid. Reis et al. [16] classified the particles of a granular 
fluid in three classes according to their shape factor. Particles with ( < ( min belong 
to class A, particles with ( min < ( < ( u belong to class B and particles with ( > ( u 
belong to class C where ( min = 1.159 and ( u = 1.25. If the number of particles in class 
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A and B are equal they observed the transition from a liquid to an intermediate phase 
(coexistence or hexatic phase could not be evaluated). Approaching the freezing point 
a sharp decline of the number of class B particles occurs in the intermediate phase. 
After crossing the freezing point the number of class B particles decreases much more 
moderately in comparison with the intermediate phase. 

Using the classification procedure described in [46J we receive an identical value for 
( min and a slightly deviating value ( u — 1.22. In accordance to the granular system the 
slope of the fraction of particle class B changes at the freezing point as shown in figure 
[TT1 On the other hand the fractions of classes A and B are equal at T 20 far away 
from the isotropic liquid - hexatic phase transition. We ascribe this different behavior 
to the structural changes which are less pronounced at the phase transitions compared 
to the hard disk or granular systems. This can been seen in figure [12] where even deep 
in the liquid state at T = 5.9 a bimodal distribution -P(C) instead of flat one exists. This 
indicates that in the liquid as well as in the crystalline state a hexagonal configuration 
of the particles is preferred. This assumption is confirmed by the fact that more than 
50% of the particles at T — 5.9 are still sixfold coordinated. For our thermal system 
with long range interaction, the phases are not identified unambiguously using shape 
factors. 
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Figure 12. The distribution of shape factors -P(C) for different interaction parameters. 
The distribution changes from a bimodal shape for T < 60.0 to a unimodal shape in 
the solid T > 60. 
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Figure 13. In 2D systems three independent Minkowski functionals exist which 
may be chosen as surface area A, the circumference U and the Euler characteristic 
X = N a - H. 



11. Minkowski functionals 

In addition to the previously described criteria we tested tentatively if any hint for 
phase transitions can be derived from the behaviour of Minkowski measures. Minkowski 
measures have been successfully adopted to specify spatial patterns, e.g. during spinodal 
decomposition [47], the evolution of galaxy clusters |48j or partial clustering in a 2D 
colloidal glass former [49] . Integral geometry offers a set of topological and geometrical 
descriptors (Minkowski functionals) to characterize spatial patterns. The operation 
of Minkowski functional V on patterns P,Q,.. have to obey three properties to be a 
morphological measure: 

(i) Motion invariance: V(gP+t) = V(P) for g and t being any rotation and translation. 

(ii) Additivity: V(P U Q) = V(P) + V(Q) - V(P n Q) 

(iii) Continuity: A slightly distortion of a pattern leads to a continuous change of V. 

According to Hadwiger [50J exist exactly D + l linear independent Minkowski functionals 
in D dimension. In case of D = 2 the Minkowski measures are related to the surface 
area A, the circumference U and the Euler characteristic \ = N a — H which is given by 
the difference of the number of connected surfaces iV a and number of holes H. 

We create a pattern for a single particle configuration by placing a cover disk 
with constant radius R at each particle position (see figure [13] a). Morphological 
information about the particle configuration is then obtained by determining the 
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Minkowski functional as a function of cover radius R. As crystals in 2D posses sixfold 
symmetry we briefly describe the behaviour of the Minkowski functionals for a perfect 
hexagonal lattice. The cover disks do not overlap for radius < R < a/2. The surface 
area and circumference are given by the area Ad(R) respectively circumference Ud{R) of 
a cover disks times the number of cover disks N and the Euler characteristic is \ = N. 
As requested by additivity for a/2 < R < v^3a/3 the surface area is given by sum of 
the area of the cover disks minus the overlapping areas which is connected with a slower 
increase of the surface area. Despite the radius is increased the circumference starts to 
decrease because the parts of the perimeters belong to the overlapping areas of the disks 
are disregarded. The Euler characteristics gets negative since the overlapping leads to 
only one connected surface while holes are formed. If the disk radius reaches a value of 
R > \/3a/3 the whole area is covered and the circumference is equal zero. The holes 
disappear, i.e. x — 1- 

In figure [13] b-d) the three Minkowski functionals b) surface area A, c) circumference 
U, and d) Euler characteristics x are shown for crystalline, hexatic and isotropic liquid 
systems. Note that A/A tota i, U/N and x/N is plotted. All Minkowski functionals 
reflect in principle the behaviour which is characteristic for the hexagonal lattice. The 
deviations from the expected curves increase with the inverse system temperature i.e. 
with thermal motion. The increase of the surface area is proportional to R 2 while the 
circumference is proportional to R and x = 1 as expected for small radii. The range 
where holes in the perfect hexagonal lattice exists are broadened to 0.35a < R < 0.8a. 
For R > 0.8a the full area is covered, U/N = and x/iV — > 0. No qualitative changes 
of the Minkowski functionals can be observed at the phase transitions. Changes of 
the Minkowski functionals which may emerge due to the emergence of defects can 
not be identified as their number density is quite small and due to strong thermal 
fluctuations of the particles. In contrast to a binary dipole-dipole system where 
structural heterogeneities have been quantified [19] Minkowski functionals are not 
applicable as order parameter to identify phase transitions temperatures for our mono- 
disperse system. 

12. Conclusion 

In a two-dimensional system of colloids with repulsive dipolar interaction several 
criteria based on structural as well as dynamical quantities were compared to identify 
phase transitions. Those criteria are the bond orientational correlation function, 
the Larson-Grier criterion, 2D dynamic Lindemann parameter, the bond-orientational 
susceptibility, the 2D Hansen- Verlet rule, the Lowen-Palberg-Simon criterion as well as 
shape factor and Minkowski functionals. A very sensitive tool to distinguish different 
symmetries ist the bond order correlation function ge(r). The transition from algebraic 
decay to exponential decay marks the hexatic to isotropic fluid phase transition. The 
bond order parameter susceptibility provides similar results for the hexatic isotropic 
transition and might be used as an alternative measure even applicable in poly- 
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crystalline samples [25] which we may name 'poly-domain hexatic' or 'polyhexaline', 
since the measure is taken above T m . 

The counterpart of g§{jc) is the density-density correlation function (eq. |5]) for the 
hexatic to crystalline transition where the crossover from quasi-long-range to short- 
range translational order marks the symmetry breaking temperature. This quantity is 
rarely used in experiment since the reciprocal lattice vector is not precise to determine 
if Mermin- Wagner fluctuations are present. Therefore we use a dynamic quantity, 
the means-squared-displacement with respect to the nearest neighbors, namely the 2D 
dynamic Lindemann parameter to identify the hexatic to crystalline transition. The 2D 
Lindemann parameter is a sensitive tool; it stays finite in the crystal but diverges in 
the fluid phase provided that the system is free of grain boundaries according KTHNY- 
theory. This should be the case for transitions with continuous character but on the 
experimental side substrate interaction, density and temperature gradients or large 
cooling rates may induce grain boundaries. 

Since the local order in 2D systems is six-fold in both, the fluid and the solid phase, 
local measures like the Larson-Grier criterion and the shape factor of voronoi cells do 
not change significantly crossing transition temperatures and are rather insensitive to 
global symmetry changes. This is at least true for our system with purely repulsive pair 
interaction where density differences do not appear in different phases. 

The Hansen- Verlet rule modified for two-dimensional systems measures the hight 
of the first peak of the structure factor. Values between S(q ) = 4.4 and S(q ) = 5.75 
are reported in simulations. In our dipolar system we determined S(qo) ~ 10 at the 
melting point. A critical value might be given for individual systems but a universal 
value should be taken with care. This is the same for the ratio of the long time versus 
short time diffusion coefficient. In 3D systems the Lowen-Palberg-Simon criterion states 
that crystallization takes place at an critical value of 0.1. In 2D values between 0.072 
and 0.099 are found in simulations, whereas in our system we found a value of 0.03. The 
discrepancies might be explained with grain boundaries where we do not know if those 
were present in the simulations. 

Finally we presented Minkowski functionals as topological measure to identify the 
phases. Whereas we found Minkowski functional to be sensitive to locally heterogeneous 
distributions of particles in a binary mixture, they were rather insensitive to global 
symmetry changes and phase transitions. 
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